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ABSTRACT 

This  thesis  discusses  the  theory  behind  the  gas  dynamic  laser. 
A  computer  program  that  solves  the  equations  is  included  along  with 
typical  outputs.   A  description  of  a  test  apparatus  is  also  presented 
with  proof  of  laser  action. 
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I.   INTRODUCTION 

The  operation  of  any  type  of  laser  system  depends  on  the  pro- 
duction of  a  population  inversion.   Descriptions  of  the  phenomena 
that  take  place  in  a  laser  may  be  found  in  [1]  and  [2].   There  are 
several  methods  of  producing  a  population  inversion,  one  of  which  is 
by  means  of  a  gas  dynamic  process.   The  gas  dynamic  method  usually 
incorporates  a  combination  of  rapid  heating  and  rapid  cooling  of  a 
gas. 

In  the  rapid  cooling  phase,  a  gas  that  is  initially  at  a  high 
temperature  and  pressure  is  exapnded  by  flow  through  a  supersonic 
nozzle,  which  causes  a  sudden  decrease  in  the  temperature  and  pressure. 
The  rapid  drop  in  temperature  is  accompanied  by  a  decrease  in  the 
population  of  the  upper  energy  levels,  each  of  which  decays  with  a 
characteristic  time  constant.   In  certain  gas  mixtures  the  time 
constants  are  such  that  a  lower  energy  level  depopulates  much  more 
rapidly  than  a  higher  energy  level.   The  result  is  a  population 
inversion  of  these  two  states.   The  method  has  been  successful  in 
producing  a  laser  output  as  shown  by  Kuehn  and  Monson  [3].   In  their 
experiment  a  large  gas  reservoir  was  heated  and  pressurized.   The  gas 
was  then  passed  through  a  nozzle  into  an  evacuated  dump  tank.   Mirrors 
were  placed  in  the  sides  of  the  nozzle  to  form  the  optical  cavity. 

It  was  the  purpose  of  this  thesis  to  build  and  test  a  gas  dynamic 
laser  utilizing  a  shock  tube  assembly  to  generate  the  initial  high 
temperature  conditions  of  reservoir  gas,  and  to  simulate  the  device  on 
the  digital  computer  according  to  the  theory  developed  by  Basov  [4,  5], 
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II.   THEORETICAL  DISCUSSION 

The  fact  that  a  gas  dynamic  laser  is  theoretically  possible  was 
established  by  Basov,  et  al.  [4].   This  work  presented  the  require- 
ments that  must  be  met  by  a  set  of  energy-level  relaxation  times  if  an 
inversion  is  to  be  established.   The  cases  of  rapid  heating  and  rapid 
cooling  were  discussed.   The  ground  states  of  a  C09-N_  mixture  were 
found  to  satisfy  the  general  requirements.   Basov,  et  al.  [5]  examined 
this  system  in  detail  and  offered  a  set  of  equations  to  describe  the 
phenomenon.   Part  of  the  results  of  their  work  is  a  system  of  equations 
that  describe  vibrational  relaxation.   These  equations  are 
dE,  r     „  E*  M„  „   ^E'  M,  2 


r        a      — 


2,0   0a  2, a/  L   2  v        2 

E0M„    x  ,E,H. 


♦«,.«..»  [-•,(t3  *•)©?)•-' 


(2) 
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Vdx~ 


r  EM  x^E'M 

=   8Z33Q3,l2[-E3C-^l+0(-9^+0 


-•Afe4  0(S)^ 


(3) 


1 "        "2 


*»HA..[-s(S+0S  +  ^Cs  +  l:Ofift) 


B  a 

dE  rxEoMox    •  x       a  /E0M.  x-| 


where 
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£-i.U-  lo"11  (.p)"1  t1/2  (i  +  i|2)  , 

£-1.095-  lo-'Ma-rtp,-1!1'^!^)"1, 

Z^1  -  1.093  ■  lo"11  (.p)-1  T1'2  (l  +  lf)\ 

_  n   10/T-l/6  /      85.08    ,    44(T\ 

Q3jl2   =  0.124T  exp^-_  +  _j, 

■,    oco    •    i«-2   m-l/6           Z'     85.08    .    440"\ 
Q3j2  =   1.858        10       T  exp  ^ T/3  +  "F"> 

Qx  2   =  0.289    •    10"6  T  exp  Qif), 

,   ,.,        1rt3  m-l/6  /'     131.055    ,    670~\ 

P2je   =6.47    ■    10     T  exp  (^-  — j^  +  — j, 

i    oo        1A4      -1/6  f      120.165    ,    612^ 

P2,a=   l'82    '    10     T  eXp  (f  —pT  +  —)' 

500    A  960    .  30 

A  =  — ,A2   =_,A3a  =  r 

Qaj3  =  0.516    .l0'7Texp(^). 
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This  gives  a  population  inversion  of 
Nr_o,  ,  -  Nr_o 


0 
where 


AN    [00  1]     [10  0]    ,    .  rz\    .  N  .        .    ,.-2  ,  .    1N-2 


cone  =  concentration  of  CO,, 


P         T 
Z  =  zr         Y  =  ^  and  e 


E 


P0        TQ      i   R0./M.  ■ 

In  the  equations,  nitrogen  is  the  "a"  gas  and  carbon  dioxide  is  the 

"3"  gas.   The  Z..  represent  the  number  of  collisions  between  a  specific 

molecule  of  type  i  and  any  j  type  molecule.   P9   is  the  probability  of 

resonant  collision  between  a  C0o  and  an  N0  molecule.   P0   is  the  corres- 

ponding  function  for  two  C0?  molecules.   Q    are  the  probabilities  of 

transition  from  one  quantum  state  to  another.   The  n  and  m  correspond 

to  the  different  energy  levels  shown  in  Figure  1.   A  listing  of  the 

(n,m)  that  goes  with  a  given  transition  may  be  found  in  [5].   The  9. 

are  the  characteristic  temperatures  of  the  corresponding  levels,  and 

hv. 
are  given  by  9.  =  — —  .   Here,  hv.  is  the  energy  of  the  given  state 

above  the  ground  state.   The  above  relations  give  the  rate  of  change 

of  energy  in  specific  energy  levels  of  C0„  and  N? .   E,  through  E 

represent  the  energies  stored  in  the  specific  energy  level.   The  initial 

9  /T 
values  of  E,  through  E  are  set  by  the  relation  E.  =  (R9./M.)/(e  i   -1) 

[5,  6]. 

A  basic  explanation  of  the  origin  of  these  equations  starts  by 

considering  the  change  of  population  of  some  state  of  one  particle  as 

the  particle  undergoes  collisions  with  another  particle.   The  rate  of 

change  of  population  would  be  given  by 
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dN.        r    N  N        -, 

-r-~-   =  Z      ~  (  S   P.   )N.  +  E   P  .  N 
dt     e,p  L  V  ,   iny   i     ,   m   nJ 
'v  n=l  n=l 

where  one  is  now  considering  an  N- level  system.   P   is  the  probability 

nm        r  .  J 

that  a  particle  will  undergo  a  transition  from  the  n   energy  level  to 
the  m   energy  level.   (Basov  uses  p   to  indicate  a  probability  of 
transition  due  to  resonant  collisions  and  Q. .  to  indicate  a  probability 
due  to  other  than  resonant  collisions.)   The  first  term  on  the  right- 
hand  side  gives  the  depopulation  of  the  i    level.   The  second  term 

introduces  the  repopulation  of  the  i    level.   Z    is  the  number  of 

e,p 

collisions  that  an  e  type  particle  will  undergo  with  a  p  type  particle. 

Note  that  only  transitions  due  to  interaction  with  another  particle 

are  being  considered  at  this  point.   Transitions  that  may  be  caused 

by  other  reasons  (spontaneous,  photon  stimulated,  etc.)  are  not  included 

These  were  assumed  negligible. 

In  order  to  describe  the  population  inversion  process,  one  must 

know  the  transition  probabilities  (p   )  and  the  collision  rates  Z 

r  'nm  e,p 

The  actual  calculations  that  must  be  carried  through  to  obtain  the  p 
will  not  be  worked  out  here.  Rather,  an  intuitive  feel  for  the  method 
used  will  be  presented.  One  starts  by  writing  Schroedinger 's  equation 
for  the  collision  process  [6].  This  will  result  in  a  wave-function 
description  of  a  particle  that  is  under  the  influence  of  an  intermole- 
cular  interaction  potential.  This  interaction  potential  is  a  function 
of  the  distance  between  the  involved  particles. 

The  resulting  wave  functions  can  be  written  down  as  some  combination 
of  the  wave  functions  of  an  isolated  particle.   The  set  of  constants 
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that  perform  this  transformation  is  directly  related  to  the  transition 

probabilities  p   .   The  relationship  between  the  elements  of  the  set 
v  rnm 

of  constants  and  the  transition  probabilities  is  known,  and  is  of 
sufficient  simplicity  as  to  allow  calculation  of  the  probabilities. 
The  origin  of  the  Z    's  is  discussed  in  [5]  where  they  are  derived 
from  the  kinetic  theory  of  rigid  molecules.   The  probabilities  given 
by  Basov  are  altered  versions  of  the  theoretical  relationships  so  that 
the  best  experimental  correlation  results. 

In  order  to  form  a  closed  system  of  equations  that  describes  the 
gas  dynamic  process,  one  must  describe  the  temperature,  pressure  and 
density  as  functions  of  time.   This  is  necessary  due  to  the  dependence 
of  the  probability  relationships  on  these  parameters.   The  quasi  one- 
dimensional  flow  equations  form  a  set  of  independent  equations  that 
describe  the  above  variables.   One  of  the  relationships  is  the  equation 
of  state : 

P  =  P  \   T  (5) 

where  R  is  the  universal  gas  constant  and  M  is  the  molecular  weight  of 
the  gas.   If  there  is  more  than  one  gas,  this  equation  can  be  separated 
into  the  partial  pressure  components.   That  is, 

1    l  M. 

l 

Assuming  steady-state  flow,  another  relationship  is  the  continuity 
equation,  which  is  given  by 

^  (pav)  =0.  (6) 

A  third  relationship  is  that  for  the  conservation  of  momentum,  as  given 
by 


—  =  -  &    .  (7) 


pV  dx     dx 
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Finally,  there  is  the  conservation  of  energy  for  adiabatic  flow.   This 
equation  takes  the  form: 

a  p 

and  in  effect  states  that  there  is  only  a  finite  amount  of  energy  that 
can  be  distributed  among  the  different  modes  of  energy  storage. 

These  eight  equations  allow  one  to  describe  the  relaxation  pehnomenon 
as  exhibited  by  an  adiabatically  cooled  gas.   A  digital  computer  program 
that  solves  the  equations  is  shown  in  the  Appendix. 

Before  describing  the  theory  of  the  gas  dynamic  process,  keeping 
the  theory  of  the  N^-CO-  gas  dynamic  laser  in  mind,  the  KL-CCL  mole- 
cular system  will  be  discussed.   The  10.6-micron  N9-C09  molecular  laser 
system  inverts  the  population  of  certain  of  the  molecular  vibrational 
energy  levels  of  the  ground  state  of  the  C0~  molecule.   The  vibrational 
energy  levels  are  split  by  the  molecular  rotational  energy.   However, 
in  the  energy-level  diagram,  for  clarity  purposes,  these  rotational 
energy  levels  are  not  shown.   The  different  vibrational  modes  of  the 
ground  state  of  the  C0„  and  N~  molecules  are  shown  in  Figure  1. 

With  reference  to  Figure  1,  E   is  the  lowest  vibrational  level  of 
the  N„  atom.   Since  N„  is  diatomic,  there  is  only  one  vibrational  degree 
of  freedom.   This  is  a  symmetrical  stretching  mode  that  can  be  visualized 
as  similar  to  two  balls  that  are  attached  by  a  spring. 

(00  1)  is  the  energy  of  the  assymetrical  stretching  mode  of  the  C0„ 
molecule.   This  is  the  case  of  the  two  oxygen  atoms  moving  in  one 
direction  while  the  carbon  atom  moves  in  the  other" direction . 

(01  0)  is  the  bending  mode  energy.   This  is  a  degenerate  state,  as 
the  molecule  can  bend  in  either  of  two  perpendicular  planes,  the  degeneracy 
being  indicated  by  the  superscript. 

17 


(10  0)  is  the  energy  of  the  symmetric  stretching  mode.   This  is 
the  same  as  in  the  case  of  N- ,  except  that  there  is  a  "stationary" 
carbon  atom  between  the  vibrating  oxygen  atoms. 

The  above  energy  levels  are  simple  harmonic  oscillator  states,  and 
therefore  the  energy  difference  between  the  nth  and  (nth+1)  energy 
levels  is  the  same,  regardless  of  the  value  of  n.   The  nomenclature 
for  the  higher  vibrational  energy  levels  is  the  same  as  above.   For 
example,  the  level  (02  0)  is  the  nondegenerate  bending  mode  that  has 
an  energy  of  approximately  twice  the  (01  0)  level.   Note  that  although 
it  is  not  shown,  the  (02  0)  would  be  right  below  the  (10  0)  level. 

There  are  other  levels  in  the  region  shown  by  Figure  1,  but  the 
levels  shown  are  those  that  enter  into  the  strong  10.6-micron  transition, 
and  are  the  most  important  ones  for  the  present  purposes. 

With  the  energy  levels  set  forth,  it  is  now  necessary  to  explain 
several  other  phenomenon  before  the  laser  can  be  completely  described. 
The  first  of  these  phenomenon  is  what  is  known  as   a  collision  of 
the  second  kind  or  a  resonant  collision.   Basically  the  resonant 
collision  is  the  interaction  of  two  atoms,  one  of  which  is  in  an 
excited  state.   The  other  atom  may  be  in  any  state.   The  result  of 
the  interaction  is  that  the  excited  atom  transfers  its  energy  to  the 
second  atom,  while  returning  itself  to  the  ground  state.   This  may  be 
written  schematically  as  [2]: 

X  +  Y*      >  X*  +  Y 
where  the  *  indicates  an  atom  in  an  excited  state.   It  has  been  shown 
that  the  probability  that  the  interaction  will  occur  decreases  as  the 
energy  gap  between  the  energy  level  of  atom  Y  and  the  closest  available 
energy  level  of  atom  X  increases.   If  the  energy  difference  is  greater 
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than  800  cm   ,  the  collision  will  effectively  not  occur  [7],   Although 
a  proof  will  not  be  presented  here,  some  insight  can  be  gained  by 
considering  the  fact  that  both  the  energy  and  the  momentum  of  the  two 
atoms  must  be  conserved,   This  criterion  is  most  readily  fulfilled  when 
there  is  little  "internal"  energy  (i.e.,  electronic,  vibrational,  or 
rotational)  that  must  be  transformed  into  kinetic  energy  of  trans- 
lation of  the  atoms. 

The  next  item  that  must  be  discussed  is  the  ftL  energy- level  diagram. 
As  seen  from  Figure  1,  for  N„  there  is  only  one  state  of  interest:   the 
lowest  vibrational  level  of  the  ground  state.   This  state  is  a  molecular 
level,  and  since  N_  has  no  net  electric  dipole  moment,  there  can  be  no 
transitions  directly  to  the  ground  state  due  to  any  electromagnetic 
interaction.   The  result  is  that  for  the  N„  molecule  to  return  to  the 
ground  state  from  the  vibrational  levels  it  must  undergo  a  collision 
with  another  particle  (other  molecules  or  the  container  walls).   This 
is  known  as  a  metastable  state,  and  is  characterized  by  a  long  lifetime. 

With  the  above  information,  the  C0?-N?  molecular  gas  dynamic  laser 
system  can  be  described.   Once  the  gas  has  been  heated,  populating  the 
upper  levels,  it  is  then  expanded  by  means  of  a  supersonic  nozzle. 
The  excited  N~  molecules  interact  with  the  C0~  molecules,  and  the  two 
types  of  molecules  undergo  resonant  collisions  which  serves  to  maintain 
the  population  of  the  00  1  level  of  the  C02 .   The  population  of  the 
other  CO^  energy  levels  depopulate  rapidly  as  both  levels  (shown  in 
Figure  1)  are  connected  to  the  ground  state  by  large  transition 
probabilities.   The  large  energy  gap  between  these  states  and  E   pro- 
hibits repopulation  due  to  resonant  collisions.   This  creates  a 
population  inversion  between  the  00°1  and  the  10°0  states.   There  are 
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other  inversions,  for  example  between  the  00  1  and  02  0  levels,  but 
the  above  inversion  is  the  strongest,  and  gives  the  10. 6-micron 
transition.   The  lasing  action  may  now  occur. 

There  is  one  other  characteristic  of  the  N„-C0?  laser  system  that 
should  be  pointed  out.   This  is  the  fact  that  any  impurities  in  the 
gases  will  cause  a  marked  decrease  in  the  output  power.   This  is  due 
to  the  resonant  transfer  of  energy  from  the  N_  molecules  to  the 
impurity  atoms,  which  robs  energy  that  would  have  otherwise  been 
pumped  into  the  inversion. 

Although  the  above  theory  does  describe  the  same  N?-C0„  laser 
system,  higher  output  powers  can  be  obtained  by  the  use  of  additives. 
In  the  case  of  a  He-ISL-CO^  laser,  a  major  factor  in  the  population  of 
the  00  1  level  appears  to  be  due  to  cascading  effects  from  higher 
energy  levels  of  both  molecules  [8,  12].   That  is,  there  is  resonant 
transfer  of  energy  between  the  C0~  and  He  at  higher  energy  levels 
(collisions  of  the  second  kind  at  upper  levels).   The  C02  molecules 
then  spontaneously  drop  into  the  00  1  level  to  set  up  the  inversion. 
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III.   EXPERIMENTAL  APPARATUS 

The  previous  section  dealt  with  the  theory  behind  the  gas  dynamic 
laser.   The  apparatus  used  to  test  the  theory  is  shown  in  Figures  2 
through  13.   The  shock-tube  facility  used  had  been  previously  built 
at  the  Naval  Postgraduate  School  and  it  is  discussed  in  detail  by 
Penaranda  [ 13 ] . 

Figure  2  shows  a  schematic  of  a  double-diaphragm  shock  tube.   The 
shock  tube  generates  a  shock  wave  by  suddenly  expanding  a  high-pressure 
gas  into  a  low-pressure  gas.   The  shock  wave  may  be  visualized  as 
formed  by  the  superposition  of  inf initessimal  compression  waves.   The 
shock  wave  is  accompanied  by  large  gradients  in  temperature  and  pres- 
sure.  The  double-diaphragm  arrangement  was  used  because  of  the  repro- 
ducability  of  run  conditions  that  it  affords.   In  the  double-diaphragm 
assembly,  one  drops  half  of  the  driver  pressure  across  each  diaphragm. 
The  diaphragms  (flat  brass  plates)  are  cut  to  withstand  this  pressure. 
However,  when  valve  #1  is  opened  the  entire  driver  pressure  is  felt 
across  diaphragm  #1  which  has  been  cut  to  rupture  well  below  this 
pressure.   Therefore,  the  shock  will  always  be  originated  by  the  same 
total  driver  pressure.   This  system  eliminates  the  necessity  that  the 
diaphragms  rupture  at  an  exact  pressure.   The  only  requirement  is  that 
they  rupture  in  the  pressure  range  between  the  full  driver  pressure 
and  one  half  that  value. 

Figure  3  shows  the  regions  of  interest  for  the  case  of  a  reflected 
shock  wave.   The  shock  wave  is  traveling  into  the  driven  gas  (region  1) 
at  a  velocity  S.   The  contact  surface  is  propagating  into  the  driven 
gas  at  a  velocity  C.   This  surface  separates  Region  2  from  Region  3. 
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Regions  1  and  2  contain  pure  driven  gas  while  Region  3  contains  pure 
driver  gas.   That  is,  there  is  no  mass  flow  across  the  contact  surface. 
The  array  of  lines  in  Region  4  represents  the  rarefaction  wave,  which 
is  produced  by  the  expansion  of  the  driver  gas.   Region  4  is  the  driver 
gas  at  rest.   Region  5  contains  pure  driven  gas  after  the  shock  has 
reflected  from  the  end  wall  of  the  shock  tube.   The  conditions  of  the 
gas  in  Region  5  were  used  as  the  initial  conditions  for  the  experiment. 

Analytical  descriptions  of  the  state  of  the  gas  in  any  of  the 
regions  may  be  obtained  by  means  of  the  conservation  of  mass,  momentum 
and  energy  equations  in  conjunction  with  the  adiabatic  gas  laws.   This 
is  done  in  (9)  and  (10).   The  results  of  interest  to  the  present  experi- 
ment are 

*5  m  V27   M2  -  (7-l)l  r<3y-l)  M2  -  2  (7-l)l 
Pl   L     <?+1>     JL   (7-l)M2  +  2      J 

h  _    f2  (7-1)  M2  +  (3-y)]  [(3-7-D  M2  -  2  (y-1)] 
Tl  [(7+1)  M]2 

where  7  is  cp/cv,  a.  is  the  speed  of  sound  in  Region  1  and  M  is  the 

mach  number  of  the  shock  wave.   The  conditions  in  Region  5  can  be 

controlled  by  variation  of  M,  P,  and  &J&-,  • 

The  remainder  of  the  physical  apparatus  is  shown  in  Figures  4 
through  12.   Figure  4  is  a  schematic  of  the  entire  setup.   The  three 
amplifiers  shown  serve  to  amplify  the  signals  from  the  thin  film 
gauges  (T.F.G.  in  the  diagram).   The  distance  between  T.F.G.  #1  and 
#2  is  5  ft.  3  in.   Therefore  by  recording  when  the  shock  passes  each 
of  the  gauges,  the  actual  shock  speed  can  be  determined.   The  schematic 
of  the  amplifiers  is  shown  in  Figure  5. 

The  biasing  circuit  for  the  detector  is  shown  in  Figure  6.   The 
detector  is  biased  at  about  370  microamps  and  gives  a  signal  of  several 
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millivolts.   The  detector  is  gold-doped  germanium  with  a  bandpass  of 
1.0  to  10.6  microns.   The  time  constant  of  the  detector  is  50 
nanoseconds. 

The  "dump  tank"  serves  as  a  reservoir  for  the  gas  after  it  has 
passed  through  the  optical  cavity  formed  by  the  mirrors.   As  will  be 
seen  in  a  latter  portion  of  this  section  ,  the   observation  time  is 
on  the  order  of  2  milliseconds.   The  particle  velocity  will  be  super- 
sonic after  the  gas  expands  through  the  nozzle.   For  a  particle 
velocity  of  4000  f ps ,  the  dump  tank  must  be  4  feet  long  to  avoid 
interference  of  the  returning  gas  with  that  gas  that  is  still  expand- 
ing through  the  nozzle.   The  gas  that  reflects  from  the  end  wall  of 
the  dump  tank  will  introduce  a  change  in  the  conditions  within  the 
gain  path.   This  change  in  the  conditions  marks  the  end  of  the  test 
time.   The  dump  tank  can  be  seen  in  Figure  7.   The  finger  like  pro- 
jection from  the  endwall  is  a  safety  release  valve,  which  is  set  to 
release  when  the  pressure  exceeds  200  psi.   This  is  well  below  the 
rated  pressures  of  all  components  in  the  system  and  will  therefore 
prevent  any  overpressures  from  building  up. 

The  mirror  and  nozzle  assembly  can  be  seen  separated  in  Figures  8 
and  9.   A  mirror  mount  is  shown  disassembled  in  Figure  10.   A  cutaway 
view  is  shown  in  Figure  11.   The  head  of  the  bolt  (shown  in  Figure  10 
with  the  mirror  attached)  is  ground  off  so  that  the  mirror  surface  is 
recessed  0.075  inch  from  the  nozzle  wall.   This  is  to  help  preserve 
the  mirrored  surface.   Original  test  runs  were  made  with  gold-coated 
glass  mirrors  (not  recessed)  and  mirror  surface  deterioration  was 
severe.   In  some  places  the  gold  was  dislodged  from  the  glass.   The 
recessing  of  the  mirrors  alleviated  this  problem. 
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The  mirror  has  a  radius  of  curvature  of  4.5  inch  (2.25-inch  focal 
length).   The  gain  path  with  the  mirrors  recessed  is  3.14  inches.   The 
mirrors  are  98%  reflective  at  10.6  microns. 

The  nozzle  is  a  linear  nozzle  with  a  half  angle  of  30  .   The  throat 
is  2  millimeters  high.   The  center  of  the  mirrors  are  1.0  inch  down  the 
nozzle  from  the  throat.   The  inside  of  the  nozzle  is  roughened  in  order 
to  minimize  internal  reflections  of  spontaneously  emitted  photons. 

The  plexiglass  shield  shown  in  Figure  12  absorbs  any  stray 
radiation.   The  end  wall  is  at  10   to  the  incident  beam  (which  enters 
from  the  left  of  Figure  12).   This  will  relfect  if   of  the  energy.   The 
reflecting  surface  is  an  IRTRAN  2  crystal.   The  hole  drilled  in  the 
plexiglass  behind  the  IRTRAN  is  to  isolate  the  IRTRAN  from  the  load 
(that  is  from  another  piece  of  plexiglass  that  will  absorb  the  trans- 
mitted beam).   The  two  metallic  bars  on  the  far  side  serve  as  counter- 
balances to  the  detector  (which  is  filled  with  liquid  N  under  operation) 
The  white  wire  carries  the  signal  while  the  black  lead  grounds  the 
detector  and  Dewar  flask. 

Figure  13(a)  shows  the  output  of  the  thin  film  gauges  of  the  shock 
tube  for  a  test  run  (0.5  m  sec/cm).   The  conditions  were  600  psi  N_ 
driving  40  mm  Hg .   ARGON  was  the  driven  gas.   The  resulting  mach  speed 
was  4.0.   An  endwall  pressure  gauge  was  substituted  for  the  nozzle  and 
its  output  is  shown  in  Figure  13(b)  (0.5  m  sec/cm).   The  flat  region  in 
the  middle  of  the  trace  is  the  observation  time  for  Region  5,  and  is 
found  to  be  on  the  order  of  2  m  sec. 
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IV.   RESULTS 

A.   THEORETICAL 

Figures  14  through  18  show  different  outputs  from  the  computer 
program.   Figure  14  shows  a  comparison  of  the  program  solutions  to 
the  predictions  by  Basov,  et  al.  [5],   The  case  considered  was  for 
an  axially  symmetric  nozzle,  with  initial  test  conditions  of  10 
atmospheres  and  1500  Kelvin.   The  two  solutions  were  found  to  differ 
by  no  more  than  3$. 

Figure  15  shows  a  comparison  of  the  linear  nozzle  as  used  in  the 
experiment  to  the  hyperbolic  nozzle  considered  by  Basov.   It  should 
be  noted  that  the  pressure  drops  faster  in  the  hyperbolic  nozzle. 
The  effect  that  this  has  on  the  population  inversion  can  be  seen  in 
Figure  16.   The  two  larger  curves  show  comparisons  of  the  computer 
predictions  to  Basov' s  predictions  for  the  hyperbolic  nozzle.   The 
smaller  population  inversion  in  the  linear  nozzle  is  caused  by  the 
higher  pressure  and  temperature  (as  compared  to  the  hyperbolic  nozzle) 
Also,  since  the  temperature  decreases  more  slowly  in  the  linear 
nozzle,  the  peak  in  the  population  inversion  occurs  farther  down  the 
nozzle  . 

Figure  17  shows  the  population  inversion  at  the  center  of  the 
optical  cavity  as  a  function  of  temperature.   There  is  some  optimum 
temperature  for  each  initial  pressure  and  C0„  concentration.   Below 
this  temperature  less  energy  is  pumped  into  the  upper  levels.   Higher 
temperatures  increase  the  gas  velocity  through  the  nozzle.   Since  the 
characteristic  relaxation  times  are  effectively  independent  of  tempera 
ture  in  this  range,  the  population  inversion  peak  occurs  farther  down 
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the  nozzle  than  the  center  of  the  optical  cavity.   This  causes  a 
smaller  inversion  at  the  center  of  the  optical  cavity. 

B.   EXPERIMENTAL 

Figure  18  shows  the  result  of  one  of  the  initial  tests.   Time 
limitations  precluded  further  experimentation.   The  setup  used  in 
this  initial  test  for  lasing  action  was  a  modification  of  that  shown 
in  Figure  4.   The  modified  setup  used  only  the  dual-beam  oscillo- 
scope.  The  oscilloscope  was  triggered  by  the  first  thin-film  gauge. 
The  second  and  third  gauges  were  fed  into  different  (i.e.,  opposite 
polarity)  differential  inputs  of  one  channel.   This  was  the  lower 
trace  in  Figure  18.   The  detector  output  was  fed  directly  into  the 
second  channel.   The  positive  jump  in  the  lower  trace  at  3.8  msec  is 
the  difference  of  the  outputs  of  the  second  and  third  thin-film 
gauges,  and  comes  just  prior  to  reflection  of  the  shock  wave  off  the 
endwall.   At  this  time  the  detector  output  shows  a  negative  step, 
which  was  caused  by  the  lasing  action.   Two  milliseconds  later  the 
conditions  of  the  gas  were  upset  and  the  observation  time  as  described 
by  the  computer  program  was  over.   This  is  noted  by  a  slight  ripple 
in  the  upper  trace  at  6.0  msec.   The  large  spikes  are  due  to  recreation 
of  favorable  conditions  by  internal  shock  reflections.   Their  large 
size  is  due  to  mixing  of  the  helium  driver  gas.   This  was  shown  by  the 
fact  that  the  spikes  disappeared  (but  not  the  original  step)  when  N_ 
was  used  as  a  driver  gas. 
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V.   CONCLUSIONS 

From  the  reproducibility  of  the  output  from  the  detector,  it  was 
concluded  that  lasing  action  was  taking  place.   However,  insufficient 
data  were  taken  to  permit  comparison  of  the  actual  laser  output  to 
the  theoretically  predicted  value.   The  only  indication  that  the  theory 
correlates  with  the  experiment  is  the  fact  that  a  population  inversion 
was  predicted  for  the  test  conditions  that  resulted  in  laser  output. 

Further  work  with  the  system,  as  shown  in  Figure  4,  is  necessary 
to  supply  the  data  required  for  full  comparisons  between  the  theory 
and  experimental  results.   Further  experimentation  with  the  laser 
should  supply  insight  into  both  the  validity  of  the  theory  and  into 
the  shock  reflection  and  gas  mixing  that  caused  the  large  output 
spikes . 
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APPENDIX 

The  eight  equations  given  in. Section  I  form  the  closed  set 
necessary  to  describe  the  gas  dynamic  laser.   The  computer  program 
at  the  end  of  this  section  solves  these  equations.   Descriptions  of 
the  variables  may  be  found  in  tables  I  and  II.   There  are  a  few 
portions  of  the  program  that  need  clarification. 

The  factor  (TT/TC)  in  statement  25  ensures  that  T  =  TT  during  the 
first  iteration.   If  T  ^  TT,  a  divergent  solution  results.   In  some 
cases  the  solution  may  converge  after  several  iterations.   The  error 
between  T  and  TT  is  small  and  is  due  to  round-off  errors  of  the 

conversion  factors  used. 

2 

The  factor  (0.98)  in  the  relation  for  CSQ  (C  )  also  ensures 

computational  stability.   That  is,  in  order  to  compute  dp/dx,  one  must 

2  2 

divide  by  (1-v  /C  ) .   At  x  =  0  (at  the  throat),  v  =  c,  and  one  would 

therefore  be  dividing  by  zero.   The  factor  0.98  ensures  that  this 
situation  does  not  occur.   The  resulting  error  is  less  than  lf>. 

The  subroutines  RKLDEQ  and  OSPLOT  may  be  found  in  the  USNPGS 
System  360  program  library.   More  information  on  them  may  be  obtained 
there. 
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TABLE    I 


Variable 

Physical  Interpretation 

Units 

Y(1),Y(2),Y(3),Y(4) 

Total  energy  in  given  level 

cm  /gr 

P5,  PI,  V(5) 

Pressure 

atmospheres 

P1AM 

Actual  driven  pressure 

millimeters 
mercury 

TT,T1,T5,T,TC 

Temperatures 

degrees  Kelvin 

Y(6),  VT 

Gas  particle  velocity 

cm/sec 

CSQ 

Speed  of  sound  squared 

2,   2 
cm  /sec 

F(1),F(2),F(3) 

,F(4) 

Population  rates 

cm  /gr/sec 

F(5) 

Pressure  rate  of  change 

atm/sec 

F(6) 

Velocity  rate  of  change 

cm/sec/sec 

PI 

Normalized  population 
inversion 

dimensionless 
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TABLE    II 


Computer  Variable 

Basov  Variable 

Y(1),Y(2),Y(3),Y(4) 

El'    V'    E3'    Ect 

Y(5) 

P 

Y(6) 

V 

F(i) 

F(i) 

T5 

T 
o 

P5 

P 
o 

CONC 

7 

hi 

hi 

Qij 

Qij 

pu 

pu 

PI 

AN/N 
o 
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FIGURE    I 
C02    AND   N2    ENERGY   LEVELS    (in   cm1) 
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FIGURE    2. 
DOUBLE     DIAPHRAGM   SHOCK  TUBE 
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Shock   Tube  At   Time  t,  >0 


Distance 


FIGURE     3 


PLOT    OF    DISTANCE     VS.     TIME    IN   SHOCK 
TUBE    AFTER  DIAPRAGM    IS  BROKEN. 
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ADJUST  NULL  BIAS  WITH  A  VTVM  BEFORE 
'  CONNECTING  OUTPUT   TO    AMPLIFIER. 


100  K  Pot  Adjusted 
For  Null  Bios 


65V.  — 


IOK 


FIGURE  6 
BIASING .  CIRCUIT 
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FIGURE    7. 
END    VIEW  OF  SHOCK  TUBE 
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FIGURE    8 
MIRROR   MOUNT 
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FIGURE    9 
NOZZLE    ASSEMBLY 
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FIGURE     10 
MIRROR  MOUNT     (DISASSEMBLED) 
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FIGURE    II. 
CUTAWAY    OF    MIRROR   MOUNT. 
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FIGURE     12 
PLEX1GLAS    SHIELD    AND    DETECTOR 
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FIGURE    13(a) 

OUTPUT    OF    T.FG.   FOR  TEST   RUN 
(0.5    m.sec./  cm) 


FIGURE    13(b) 

OUTPUT   OF  ENDWALL  PRESSURE  GAUGE 
(0.5  m  sec ./ cm) 
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Detector 


T.F.G. 


2  msec/cm 
Lower  and  upper  trace:  10  mv/cm 

Run  conditions: 

Driven  gas:  200  ram  (Hg)  of  20,7%   CO  in  W2 

Driver  gas:  600  psi  of  40$  Argon  in  Helium 


FIGURE  18 
DETECTOR  OUTPUT  SHOWING  LASER  ACTION 
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